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O 

' ning of the precision era of Cosmic Microwave Background anisotropy (CMB) measurements. We 

investigate the observational constraints from the current CMB anisotropy measurements on the 
simplest inflation models, characterized by a single scalar field </!>, in the parameter space consist- 
ing of scalar spectral index ns and tensor/scalar ratio r. If we include constraints on the baryon 
density from big bang nucleosynthesis (BBN), we show that the favored inflationary models have 
negligible tensor amplitude and a "red" tilt, with a best flt of ns — 0.93, which is consistent with 
the simplest "small-fleld" inflation models, but rules out large-field models at the la level. Without 
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1/^ ' including BBN constraints, a broader range of models are consistent with the data. The best flt 

m 

' and hybrid scenarios. Large-field models (such as those typical of the chaotic infiation scenario) 

■ with tilt ns < 0.9 are strongly disfavored in all cases. 

O ■ 98.80.Cq,98.70.Vc,98.80.Es; SNS-PH-00-12 

I. INTRODUCTION 

It is commonly believed that the Universe underwent an early era of cosmological inflation. The flatness and 
the horizon problems of the standard big bang cosmology are elegantly solved if, during the evolution of the early 
Universe, the energy density is dominated by some vacuum energy and comoving scales grow quasi-exponentially. The 
prediction of the simplest models of inflation is a flat Universe, i.e. ^tot — 1 with great precision. 

Inflation ^ has also become the dominant paradigm for understanding the initial conditions for structure formation 
and for Cosmic Microwave Background (CMB) anisotropy. In the inflationary picture, primordial density and gravity- 
wave fluctuations are created from quantum fluctuations "redshifted" out of the horizon during an early period of 
superluminal expansion of the universe, where they are "frozen" as perturbations in the background metric [^-^. 
Metric perturbations at the surface of last scattering are observable as temperature anisotropy in the CMB. 
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The first and most impressive confirmation of the inflationary paradigm came when the CMB anisotropics were 
firmly detected by the COBE satellite in 1992 Subsequently, it became clear that the measurements of the spec- 

trum of the CMB anisotropy can provide very detailed information about the fundamental cosmological parameters 
I pof and other crucial parameters for particle physics This era of precision CMB anisotropy measurements has 
just begun with the results of two balloon-borne experiments, Boomerang ]l^ , p^ ~ which was flown over Antarctica 
in 1999 - and MAXIMA-1 ||l^,|l^. The observed first acoustic peak at ^ ~ 200 (already detected by previous experi- 
ments, see e.g. |l^Jl^) indicates that the curvature of the Universe is consistent with flatness, thus providing another 
convincing confirmation of the fact that the early Universe experienced a period of superluminal acceleration. 

Despite the simplicity of the inflationary paradigm, the number of inflation models that have been proposed in 
the literature is enormous ||l^. Models have been invented which predict non-Gaussian density fluctuations 
isocurvature fluctuation modes ||2^, and cosmic strings At present, inflation plus topological defects is consistent 
with the Boomerang data and the evidence for non-Gaussianity in the COBE maps is tantalizing but not 
compelling. Furthermore, a small class of models predicting f^o < 1 (open inflation models) have not survived the 
current data. 

The goal of this paper is to discuss the capability of the existing CMB anisotropy data to discriminate among 
the very simplest type of inflationary models: inflation involving a single scalar fleld (the inflaton) producing a flat 
Universe consistent with the observed first acoustic peak at £ ~ 200. The single-field models are certainly the most 
popular among models for cosmological inflation because they are flrmly rooted in modern particle theory and have 
usually supersymmetry as a crucial ingredient [ p^ . 

For single-field inflation models, the relevant parameter space for distinguishing among models is the plane defined 
by the scalar spectral index ns and the ratio of tensor to scalar fluctuations r. Recent work has shown that there is a 
kincmatically favored region in the r — ns plane, independent of the details of the model |2j]. Within that "attractor" 
region, different models make distinct predictions for r and ns- In this paper we show that the Boomerang and 
MAXIMA-1 data sets allow us to significantly constrain the allowed region in the r — ns plane and begin to rule out 
particular models of infiation. Such constraints are sensitive to our assumptions about other cosmological parameters, 
in particular the baryon density f2b and the reionization optical depth Tc. If we allow both of these parameters to 
vary freely, a fairly large region in the r — ns plane is consistent with observation. A scale-invariant spectral index 
715 = 1 is favored, and the models which predict both a tilted spectrum 715 < 1 and large tensor modes r > 0.5 
are strongly disfavored. If, conversely, we choose to constrain ilb using constraints from primordial nucleosynthesis, 
the allowed region in the parameter plane is tightened considerably. If we additionally assume negligible reionization 
optical depth, the data significantly favor a tilted spectrum ns < I and negligible tensor modes. 

The paper is organized as follows. In Section II we describe the classes of infiation models we consider and their 
predictions for cosmological parameters which can be distinguished via CMB observations. In Section III we discuss 
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the methods used for the CMB analysis. In Section IV we present the constraints from the combined Boomerang and 
MAXIMA-1 data sets. In Section V we discuss some general conclusions. 



II. THE INFLATIONARY MODEL SPACE 



In single-field models of inflation, the inflaton "field" need not be a fundamental field at all, of course. Also, some 
"single-field" models require auxiliary fields. Hybrid inflation models |25-2^], for example, require a second field to end 



inflation. What is significant is that the inflationary epoch be described by a single dynamical order parameter, the 
inflaton flcld. The metric perturbations created during inflation are of two types: scalar, or curvature perturbations, 
which couple to the stress-energy of matter in the universe and form the "seeds" for structure formation, and tensor, 
or gravitational wave perturbations. Both scalar and tensor perturbations contribute to CMB anisotropy. Scalar 
fluctuations can also be interpreted as fluctuations in the density of the matter in the universe. Most (but not all) 
inflationary models predict that these fluctuations are generated with approximately power law spectra 



Ps (k) oc fc"«-\ 

PT{k)(xk"^. (1) 

The spectral indices ns and tit are assumed to vary slowly or not at all with scale: 

~ (2) 

In actuality, for a spectral index sufficiently higher than unity, the power-law approximation becomes less and less 
accurate [|4 28 1 . Some particular inflation models predict running of the spectral index with scale ||29|-|3^ or even sharp 



features in the power spectrum |37|. The observational consequences of scale dependence of the spectral index are 
considered in Ref. |]3l[| . We do not consider such models here. In principle, then, we have four independent observables 
with which to test inflation. One can be removed by considering the overall amplitude of CMB fluctuations: if the 
contribution of tensor modes to the CMB anisotropy can be neglected, normalization to the COBE four-year data 
gives p8|,^ Pg^^ = 4.8 x 10^^. We then have three relevant observables: the tensor scalar ratio deflned as a ratio of 
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quadrupole moments 
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r ^ --^ (3) 



^Scalar 

and the two spectral indices ns and nr- 

Even restricting ourselves to a simple single-field infiation scenario, the number of models available to choose from 
is large It is convenient to define a general classification scheme, or "zoology" for models of inflation. We 

divide models into three general types: large-field, small-field, and hybrid, with a fourth classification, linear models, 
serving as a boundary between large- and small-field. A generic single-field potential can be characterized by two 
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independent mass scales: a "height" A^, corresponding to the vacuum energy density during inflation, and a "width" 
/i, corresponding to the change in the field value /S.4> during inflation: 



V{4>)^^'f{j^- (4) 



Different models have different forms for the function /. The height A is fixed by normalization, so the only free 
parameter is the width /i. To create the observed flatness and homogeneity of the universe, we require many e-folds 
of inflation, typically N ~ 50. This figure varies somewhat with the details of the model. During inflation, scales 
smaller than the horizon are "redshifted" to scales larger than the horizon. A comoving scale k crosses the horizon 
during inflation N (k) e-folds from the end of inflation, where N (k) is given by 

Here Vk is the potential when the mode leaves the horizon, Ve is the potential at the end of inflation, and prh is 
the energy density after reheating. Scales of order the current horizon size exited the horizon at TV (fc) ~ 50 — 70. 
In keeping with the goal of discussing the most generic possible case, we will allow N to vary within the range 
50 < A^ < 70 for any given model. Useful parameters for our purpose are the so-called slow roll parameters e and ry, 
which are related to the potential as |M : 



Mil (V'{<j,)^^ 



and 



Stt 



167r \ V{(I)) 



V{<P) 2\V{cj>) 



(6) 



(7) 



We can write the observable parameters r, ns and nT in terms of the model parameters e and 77 ||4l|| : 

r = 13. 6e, 
= 1 - 4e + 277, 

TIT ^ -2e. (8) 
Note that r, ns and ut are not independent. The tensor spectral index and the tensor/scalar ratio are related as 

nr = -^r, (9) 
known as the consistency relation for inflation. (This relation holds only for single-field inflation, and weakens to an 



inequality for inflation involving multiple degrees of freedom |42-44|.) The relevant parameter space for distinguishing 
between inflation models is then the r — ng plane. Different classes of models are distinguished by the value of the 
second derivative of the potential, or, equivalently, by the relationship between the values of the slow-roll parameters 
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e and r^. Each class of models has a different relationship between r and ns- For a more detailed discussion of these 
relations, the reader is referred to Refs. [ ^7|j4^ . It should be noted that the accuracy of the slow roll approximation 
itself begins to break down for a spectral index significantly far from ns — 1 p8| . First order in e and r] is sufficiently 
accurate for the purpose of this analysis. 

A. Large-field models: — e < 77 < e 

Large-field models are potentials typical of the "chaotic" inflation scenario , in which the scalar field is displaced 
from the minimum of the potential by an amount usually of order the Planck mass. Such models are characterized 
by V" {(p) > 0, and —€< rj < e. The generic large-field potentials we consider are polynomial potentials V {(p) — 
A'*(0//i)^, and exponential potentials, V {(p) = A^exp For the case of an exponential potential, ¥{(/)) oc 

exp the tensor/scalar ratio r is simply related to the spectral index as 

r = 7(l-ns). (10) 

This result is often incorrectly generalized to all slow-roll models, but is in fact characteristic only of power-law 
inflation. For inflation with a polynomial potential V {(p) (x (p'^ , we again have r oc 1 — ns, 

r^7(^){l-ns). (11) 



so that tensor modes are large for significantly tilted spectra. 

B. Small-field models: 77 < —e 

Small-field models are the type of potentials that arise naturally from spontaneous symmetry breaking (such as the 
original models of "new" inflation [^,0 ) and from pseudo Nambu-Goldstone modes (natural inflation ) . The field 
starts from near an unstable equilibrium (taken to be at the origin) and rolls down the potential to a stable minimum. 
Small-field models are characterized by V" {(p) < and rj < — e. Typically e (and hence the tensor amplitude) is close 
to zero in small-field models. The generic small-field potentials we consider are of the form V {(p) = [1 — {(p/ fjif], 
which can be viewed as a lowest-order Taylor expansion of an arbitrary potential about the origin. The cases p — 2 
and p > 2 have very different behavior. For p = 2, 

r = 7(l-ns)exp[-l-Ar(l-ns)], (12) 
where iV is the number of e-folds of inflation. For p > 2, the scalar spectral index is 



'^The designations "small-field" and "large-field" can sometimes be misleading. For instance, both the 7?^ model [glj and the 
"dual inflation" model are characterized by A<j!> ~ mpi, but are "small-field" in the sense that 77 < < e, with ns < 1 and 
negligible tensor modes. 
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N \p-2 

independent of r. Assuming jj, < Mpi results in an upper bound on r of 



r^s^l-i-A^]^ (13) 



""^ N{p-2) [np{p-2)) ■ ^ ' 



C. Hybrid models: < e < r; 

The hybrid scenario j2^-|2^ frequently appears in models which incorporate inflation into supersymmetry. In 
a typical hybrid inflation model, the scalar field responsible for inflation evolves toward a minimum with nonzero 
vacuum energy. The end of inflation arises as a result of instability in a second field. Such models are characterized 
by V" {(/)) > and < e < 77. We consider generic potentials for hybrid inflation of the form V {(p) = A'' [1 + {(f)/ p)^] . 
The field value at the end of inflation is determined by some other physics, so there is a second free parameter 
characterizing the models. Because of this extra freedom, hybrid models fill a broad region in the r — ns plane (Fig. 
|l|). There is, however, no overlap in the r — ns plane between hybrid inflation and other models. The distinguishing 
feature of many hybrid models is a blue scalar spectral index, 715 > 1. This corresponds to the case 77 > 2e. Hybrid 
models can also in principle have a red spectrum, n < 1. 

D. Linear models: — —e 

Linear models, V (</>) oc (p, live on the boundary between large-field and small-field models, with V" (</>) = and 
r] = — e. The spectral index and tensor/scalar ratio are related as: 

r = ^(l-ns). (15) 

This enumeration of models is certainly not exhaustive. There are a number of single-field models that do not 
fit well into this scheme, for example logarithmic potentials V (cj)) oc In (0) typical of supersymmetry p8[ |. Another 
example is potentials with negative powers of the scalar field V {(f) oc used in intermediate infiation ||53|] and 
dynamical supersymmetric inflation p2| , ^ . Both of these cases require and auxilliary field to end infiation and are 
more properly categorized as hybrid models, but fall into the small-field region of the r — ns plane. However, the 
three classes categorized by the relationship between the slow-roll parameters as — e < 77 < e (large-field), rj < —e 
(small-field, linear), and < e < 77 (hybrid), cover the entire r — ns plane and are in that sense complete.^ Figure ^ 
J47[ shows the r ~ ns plane divided up into regions representing the large field, small-field and hybrid cases. 



■^Ref. [|48| incorrectly specified < 77 < e for large-field and 77 < for small-field. 
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III. CMB ANALYSIS 



In this section we discuss the abihty of current CMB data to place constraints on the inflationary parameter space. 
In particular, we wish to determine which regions of the inflationary model space are ruled out by the Boomerang and 
MAXIMA- 1 data, and which regions are consistent with the observed data. We perform a likelihood analysis over 
multiple cosmological parameters and project likelihood contours onto the r ~ ns plane. The choice of parameters to 
be varied in the analysis is crucial, and different assumptions result in different constraints on the inflationary model 
space. In addition to ns and r, we also vary the Hubble constant /i, the reionization optical depth Tc, the amplitude 
of fluctuations, Cio, (in units of C^^^^), the cold dark matter density f^M, the vacuum energy density J^a, and the 
baryon density JIm, subject to the constraint of a flat universe fitotai f^B + f^CDM + f^A = 1, consistent with the 
prediction of inflation. We also assume the consistency relation between tlt and r. We fix the number of neutrino 
species at three and we take a gaussian prior on the Hubble constant: h = 0.65 ± 0.2. Note that this approach is 
different from many CMB parameter estimation programs, which attempt to constrain cosmological parameters in a 
model independent fashion. Since our object here is to falsify models, we include as many model dependent constraints 
as possible and examine the consistency of the data with the predictions of the models. For example, this allows us 
to avoid questions of parameter mis-estimation arising from features in the primordial power spectrum js^ . 

The Boomerang and Maxima power spectra are estimated in 12 and 10 bins respectively, spanning the range 
2^ < t < 785. In each bin, the spectrum is assigned a flat shape, 1)Ci/2'k = Cb- Following j5^ we use the offset 
lognormal approximation to the likelihood L. In particular we define ^: 

- 21nL = {D%' - D'^i)MBB'{D'^, - D%-,), (16) 
=lniC§ +xb), (17) 

Mbb' = {C%^ + xb)Fbb'{C^b' + XB'), (18) 

where (C'j^) is the theoretical (experimental) band power, xb is the offset correction and Fbb' is the Gaussian 
curvature of the likelihood matrix at the peak. Proceeding as in Refs. we assign a likelihood to each point in 

the parameter space and we find constraints on all the parameters by finding the remaining "nuisance" parameters 
which maximize it. 

Usually there are two maximization procedures. One is based on a search algorithm through the second derivative of 



the likelihood matrix |57|. In this approach the C^'s are computed on the way, without sampling the whole parameter 



^The cross-correlation between bandpowers and the offset lognormal corrections for the Boomerang and Maxima experiments 
are not yet public available. We tested the stability of our result including a 10% correlation and with different lognormal 
distributions. We obtained near identical results with the analysis of where the comparison was possible (f2 = 1, r = 0). 
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space. The second approach is based on building a database of Ci's on a discretized grid of the parameter space. 
The maxima values are then obtained from the likelihood computed on the grid [ p8|j59| ]. We adopted the database 
approach mainly for its robustness and the possibility to use a 'true' marginalization method'. The bulk of our results 
were evaluated on a database of models sampled as follows: 

= ^CDM + f^B = (0.10, 0.15, 1.1) 
Ob = (0.015,0.030, ...,0.20) 
VLtv (0.0,..., 0.95) 
n, = (0.50, 0.52,..., 1.48, 1.50) 

(0.25,0.30, ...,0.95) 
Tc = (0.0,0.2,0.4) 
fit — — r/6.8 
r continuous 
Cio continuous. 

Two parameters in particular require careful consideration: the physical baryon density Q.\yh? and the reionization 
optical depth Tc. Our constraints on the inflationary parameter space are strongly dependent on our choice of priors 
for Q,hh? ■ In particular, the suppressed second peak seen in the data requires either a baryon density too high to be 
consistent with big-bang nucleosynthesis (BBN) constraints or a significant "red" tilt, ns < 1, if standard BBN is 
assumed. For completeness, we will consider separately 3 cases: no prior constraints on fl^h^, 0.005 < n^,h^ < 0.014 
and 0.016 < il\^h^ < 0.021 as suggested by the dichotomy in high and low deuterium measurements from distant 
quasars absorption line systems (see, e.g. Rcf. p0| , |6l[ ). The second important parameter is the reionization optical 
depth Tc- It is known that there is a strong degeneracy between Tc and the tensor/scalar ratio r, so that allowing 
for significant reionization strongly degrades the ability to distinguish tensor modes ji^ . Theoretical predictions, 
based on the Press-Schecter formalism, show that reionization most probably occurred after redshift z ^ 50 ||6^, with 
z ^ 10 — 20 for A CDM flat models []63| , |64[ . This roughly corresponds to an optical depth < Tc < 0.2 |^ . Again, 
we consider separately 3 cases: no reionization, Tc — 0.2 and Tc — 0.4. 

IV. RESULTS 

In the r — ng-plane, the parameter space is divided into regions for small-field, large-field and hybrid models as in 
Figure 1. Figure 2 shows the la {6x^ — 2.3), 2a {6x^ = 6.0), and 3a {6x^ — 9.2) contours on the r — ns plane for the 
case of no reionization and no BBN prior. The various inflation models described in the last section and in Figure 1 
are now plotted as labeled lines on the graph. (Note that the entire region to the right of the case of the exponential 
potential is consistent with hybrid inflation models). The lines are obtained by varying a given parameter within a 
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class of models, e.g. for large-field models with a power-law potential cj)'' what varies is the power p. 

The best model is nearly scale invariant, and the tensor/scalar ratio r is only weakly constrained. Marginalizing over 
r [us) we obtain also the 1 — ct constraint : 0.94 KUg < 1.07 (r < 0.65). Hybrid and small-field models are consistent 
with the data, but strongly tilted large-field models are in conflict with the CMB constraint to high significance. The 
3cr contour on the V {(f)) oc 0^ models result in a constraint that p < 5 for = 50 and p < 8 for N ~ 70. In the 
power-law inflation case, V {(f)) oc exp((/)//i), the 3cr constraint corresponds to a lower limit /i > 0.75Afpi. The la 
contour is marginally consistent with p — 2 for = 50 and requires p < 3 for = 70 in the case of a polynomial 
potential, and constrains fi > Mpi in the power-law inflation case. 

Figure 3 shows a similar plot when BBN constraints are included, still assuming negligible reionization. The 
contours in r are significantly tighter than in the case with no BBN prior, and the best fit model is shifted toward 
tilted models: 0.88 < ns < 0.98 and r < 0.17 for the low deuterium (LD) case ; 0.87 < ns < 0.97 and r < 0.16 for the 
high deuterium (HD) case, with the scale invariant case just consistent with the data at 2a. (Fig. 3 shows the high 
deuterium case.) To la, large-field models are excluded entirely. The only models within la are linear and small-field 
inflation. To 3a, the case of a polynomial potential is constrained to p < 6 for = 50 and p < 8 for A^ = 70. 

Figure 4 shows the constraints when a strong constraint on fib from BBN is assumed (left to right) and when 
the optical depth Tc is increased. Fixing the optical depth to Tc ^ 0.2 shifts the likelihood towards "blue" tilted 
models {ns > 1) making the scalar invariant models consistent with the data even when the BBN constraint are 
assumed (1.02 < ns < 1.28 for no constraint, 0.95 < ns < 1.05 for both BBN constraints, all at 68% c.l.). Increasing 
this parameter up to Tc — 0.4 moves the 68% likelihood contours in the hybrid models region (1.15 < ns < 1.30, 
1.06 < ns < 1.14 for LD and 1.00 < ns < 1.08 for HD). These models arc disfavored by the dataset itself, being 
the best flt model at Ax^ ^ 3 from the corresponding best fit model with Tc ~ 0. Note in particular that large-field 
models with tilt n^ < 0.9 are strongly disfavored in all cases. Removing the consistency relation and fixing nt = 
does'nt affect this conclusion and has small effect on the overall result: for no BBN constraint and Tc = we found 
r < 0.67 and 0.96 < < 1.14 (marginalized over r). 

V. CONCLUSIONS 

With the release of data from the Boomerang and MAXIMA-1 balloon fiights, it is possible for the first time to place 
significant constraints from the cosmic microwave background on the space of possible inflation models. The data in 
the region £ > 300, in particular, break the degeneracy between ns and r on the first peak and produce a more stringent 
upper limit on r for hybrid models respect to previous analysis ]6^-|6^. In this paper we consider the observational 
constraints on the "zoo" of simple infiation models characterized by a single scalar field (/) in the parameter space 
consisting of scalar spectral index ns and tensor/scalar ratio r. Different classes of models make distinct predictions 
for the relationship between ns and r and can in principle be differentiated by use of CMB observations. While the 
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current data is not good enough to rule out entire classes of models, it is sufficient to place significant constraints 
on model parameters. Combined with constraints on the baryon density from Big-Bang nucleosynthesis, we conclude 
that the favored inflationary models have negligible tensor amplitude and a "red" tilt, with a best fit of r < 0.17 
and ns = 0.93 ±0.05 (assuming low deuterium and negligible reionization) , consistent with the simplest "small-field" 
inflation models such as those arising from spontaneous symmetry breaking |^,^ or from pseudo Nambu-Goldstone 
bosons (natural inflation |^2|). Models with strong reionization favor higher values of the spectral index consistent 
with hybrid inflation models. Without including BBN constraints, a broad range of models are consistent with the 
data, including large-field and hybrid scenarios, and the favored model is nearly scale-invariant, ns — I. Future 
observations promise to significantly narrow the allowed regions of the parameter space, and will potentially make it 
possible to rule out entire classes of infiation models. 
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FIG. 1. The parameter space divided into regions for small-field, large-field and hybrid models. The linear case is the dividing 
line between large- and small-field. 
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FIG. 2. CMB constraints and inflation models for Tc = and no BBN prior. The aUowed contours are quite large but still 
exclude a significant portion of the inflationary model space. 
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